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Semianalytical  expressions  are  developed  for  creeping  flame  spread  parallel  to  a  fuel  cylinder  axis  so 
that  spread  rates  measured  or  predicted  for  flat  surfaces  over  the  same  material  can  be  directly  used  for 
cylindrical  surfaces.  Two  phenomena  are  identified  that  cause  an  increased  flame  spread  rate  in  the  cylin¬ 
drical  over  the  flat  surface  geometry:  (1)  increased  heat  transfer  from  the  gas  phase  and  (2)  faster  tem¬ 
perature  increase  of  the  solid  phase.  Analytical  expressions  from  approximate  solution  to  a  conserved  scalar 
equation  are  developed  for  the  flame  spread  rate  by  simplifying  the  flowfield  using  an  Oseen  approximation 
for  the  opposing  flow.  The  spread  rate  expression  for  the  fuel  cylinder  is 

pMT,  -  TJVfL,^  (^1  -  =  C(Lg/RJ£FL/ln(l  +  C  -  L^a{Tt  -  if} 

in  which  the  terms  containing  the  constant  C  account  for  the  enhanced  gas-to-surface  heat  transfer  because 
of  the  cylindrical  curvature,  and  those  containing  L,y,  the  heated  layer  depth  in  the  solid,  account  for  a 
reduction  in  the  solid  volume  preheated  in  the  cylindrical  compared  to  the  flat  geometry.  The  expression 
is  tuned  by  comparison  with  complete  numerical  solutions  to  the  flame  spread  problem  from  which  the 
flame  energy  Ejtl  is  determined  from  the  flat  surface  geometry  and  the  constant  C  chosen  from  heat 
transfer  correlations.  Results  compare  favorably  with  numerical  solutions  for  cylindrical  spread  in  forced 
and  natural  flows  and  microgravity  and  with  experiments  on  dovraward  flame  spread  on  cylindrical  rods 
in  normal  gravity  and  microgravity. 


Introduction 

Opposed-flow  flame  spread  over  flat  surfaces,  also 
called  creeping  flame  spread,  has  been  studied  ex¬ 
tensively  during  the  past  three  decades  [1-3].  Al¬ 
though  the  influence  of  parameters  such  as  ambient 
oxygen  content,  pressure,  and  fuel  thickness  is  well 
understood,  few  modeling  studies  exist  on  the  effects 
of  fuel  geometry  (e.g.,  the  cylindrical  geometry  con¬ 
sidered  here).  Ref.  [4]  is  one  exception  in  which  ra¬ 
dial  convergence  of  the  conduction  heat  flux  from 
the  gas  to  the  solid  was  considered,  but  radiation. 


important  at  microgravity  and  included  here,  was  ig¬ 
nored.  Fire  spread  in  a  cylindrical  geometry,  such  as 
along  cable  insulation,  is  a  serious  fire  safety  concern 
in  microgravity  [5] .  Previous  investigations  and  mea¬ 
surements  of  creeping  flame  spread  along  cylindrical 
rods  include  the  work  of  Sibulkin  and  Lee  [6]  and 
Fernandez-Pello  and  Santoro  [7]  at  normal  gravity 
and  experiments  in  low  gravity  in  the  European 
space  program  [8]. 

The  present  analytical  and  numerical  modeling 
was  initiated  in  support  of  experiments  on  flame 
spread  along  cylindrical  rods  or  tubes  as  part  of  the 
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FIRE — Ignition  and  Flame  Spread 


Fig.  I.  Schematic  of  the  flame 
spread  problem. 


reflight  of  the  Solid  Surface  Combustion  Experi¬ 
ment  (SSCE)  [9-11].  The  objectives  here  are  to  (1) 
present  the  physics  for  creeping  flame  spread  along 
a  cylindrical  geometry,  (2)  develop  semianalytical  ex¬ 
pressions  supported  by  numerical  modeling  that  re¬ 
late  flame  spread  in  the  cylindrical  geometry  with 
flame  spread  in  the  flat  geometry,  and  (3)  compare 
predictions  with  existing  data  on  cylindrical  surfaces 
in  normal  gravity  and  microgravity  [6,7,12]. 


Background 

Review  of  the  Flat  Surface  Energy  Balance 

Consider  the  control  volume  near  the  flame  front 
shown  in  Fig.  1  in  flame-fixed  coordinates  with  the 
oxidizer  approaching  the  flame  with  speed  =  Vg 
+  Vf  and  the  fuel  approaching  the  flame  at  Vf,  where 
Vg  is  the  opposing  flow  velocity,  and  Vf  is  the  spread 
rate.  In  the  following  discussion,  the  opposing  flow 
velocity  profile  is  flat  (Oseen  approximation),  and 
gaseous  thermal  properties  are  constant,  indepen¬ 
dent  of  temperature.  A  methodology  to  account  for 
variable  properties  is  outlined  later.  The  solid  is  ther¬ 
mally  thick,  and  the  thermally  thin  situation  is  ob¬ 
tained  as  a  limiting  case. 

The  heat  transfer  rate  per  unit  sample  width  from 
the  flame  near  the  leading  edge  expended  in  the 
sensible  heating  of  the  solid  phase  from  ambient  to 
the  vaporization  temperature  (i.e.,  T„  to  Tf  as  it  ap¬ 
proaches  the  flame  is 

ACs(Tv  —  TfVfh^  =  £fl  “  l^RAD  (1) 

The  length  scale  in  the  solid  phase  is  the  normal 
depth  over  which  the  solid  is  preheated  as  it  ap¬ 
proaches  the  flame  front  where  the  flame  heat  flux 
extends  over  the  gas-phase  thermal  length  Lg,  which 
equals  Og/V^,  with  Og  the  thermal  diffusivity  of  the 
gas,  and  and  the  solid  density  and  specific  heat, 
respectively.  equals  [9] 

Jajt  =  Ja^L^/Vf 

whereas  the  axial  length  scale  in  the  solid,  is 


aJVf,  with  t  the  time  the  solid  takes  to  pass  through 
the  flame  region  of  length  Lg.  For  L^y  >  t,  the  fuel 
thickness,  the  fuel  is  thin,  and  L^y  is  replaced  by  t. 

The  right-hand  side  of  equation  1  contains  the  en¬ 
ergy  conveyed  by  the  flame  to  the  solid,  LpL,  includ¬ 
ing  convection  and  gas-phase  radiation,  for  increas¬ 
ing  the  sensible  enthalpy  of  the  solid,  reduced  by  the 
energy  lost  through  surface  reradiation,  £rad>  and  it 
takes  the  form  [13] 

£fl  =  -  KB  -  r)L]Lg  =  i 

Cp  E^p 

[(B  -  r)L]Lg  =  ^  (B  -  r)L  (2) 
Cp 

where  the  mass  transfer  number 

B  =  [(/A  -  /r)Yo^o  -  Cp  (T,  -  T„)]/L 
with  the  heat  of  combustion,  Y„  the  oxygen  mass 
fraction,  L  the  heat  of  vaporization,  Cp  the  gas  spe¬ 
cific  heat,  r  the  stoichiometric  fuel/ oxidizer  ratio,  kg 
the  gas  thermal  conductivity,  the  convective  heat 
transfer  coefficient,  Xa  the  combustion  efficiency, 
and  Xr  the  radiative  loss  fraction.  Erad  is  equal  to 
the  rate  of  heat  transfer  per  unit  width  over  a  gase¬ 
ous  thermal  length  that  also  determines  the  region 
along  the  solid  heated  by  the  flame  leading  edge, 
that  is,  Erad  =  frrLg  =  with  e  the 

surface  emittance,  set  equal  to  unity  here,  and  a  the 
Stefan-Boltzmann  constant. 

The  energy  balance  of  equation  1,  equivalent  to 
the  de  Ris  expression  except  for  a  constant  [14],  to¬ 
gether  with  the  definitions  for  L„,  L  Epp,  and  Erad 
allow  the  flame  spread  rate  to  be  determined.  The 
formulation  can  also  be  generalized  for  different 
profiles  of  opposing  flow  by  noting  that  the  flame 
spread  rate  can  be  interpreted  to  depend  on  two 
parameters:  the  flame  energy,  Epp,  and  the  gaseous 
thermal  length.  Eg. 

There  are,  for  all  intents  and  purposes,  three 
spread  rate  regimes:  (1)  a  thermal  regime  at  mod¬ 
erate  opposing  flow  velocity  in  which  radiation  is  un¬ 
important,  and  the  flame  energy  becomes  a  property 
of  the  material  and  its  environment,  independent  of 
Eg  for  a  thin  fuel;  (2)  a  high  opposing  flow  velocity 
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regime  in  which  radiation  is  unimportant,  but  EpL 
depends  on  a  Damkohler  number  and  hence  Lg, 
even  for  a  thin  fuel,  and  the  shape  of  the  velocity 
profile;  and  (3)  a  low  opposing  flow  velocity  regime 
in  which  EpL  depends  on  Eg,  at  least  to  an  extent, 
through  the  dependence  of  the  radiation  on  Eg.  As 
Vg  — >  0  in  the  quiescent  microgravity  environment 
(i.e.,  the  limiting  case  of  the  low  opposing  flow  ve¬ 
locity  regime),  Vf  is  the  only  velocity  that  appears, 
and  EpL  depends  primarily  on  material  and  environ¬ 
mental  properties.  With  the  exception  of  the  high 
opposing  flow  velocity  regime,  we  can  view  EpL  as  a 
material  property  that  may  be  determined  by  mea¬ 
surement  [13]  or  from  numerically  determined 
spread  rates  in  conjunction  with  equation  1  account¬ 
ing  for  variable  properties  and  the  opposing  flow  ve¬ 
locity  profile,  which  contributes  to  determining  h^. 

Effects  of  Cijlindrical  Geometry 

Application  of  equation  1  can  be  extended  to 
creeping  flame  spread  along  a  cylinder  by  emphasiz¬ 
ing  the  similarities  and  differences  between  the  flat 
and  cylindrical  geometry  with  the  same  opposing 
flow  velocity: 

1.  The  gaseous  length  scale.  Eg,  remains  the  same. 

2.  The  flame  energy  transferred  to  the  solid  in¬ 
creases  because  the  heat  transfer  coefficient,  h^, 
increases  with  the  curvature  [4,15]. 

3.  The  sensible  heat  required  to  raise  the  solid  sur¬ 
face  temperature  from  ambient  to  the  pyrolysis 
value  decreases  compared  to  a  flat  surface  of  the 
same  thickness  as  the  radius  because  the  cylin¬ 
drical  surface  seen  by  the  flame  is  convex. 


Model  Development 


Flame  Energy  in  the  Cylindrical  Geometry 

It  is  well  known  that  heat  fluxes  in  laminar  flow 
along  a  cylinder  increase  as  the  radius  decreases  (see 
Ref.  [15]  for  the  hot  cylinder  and  Ref.  [16]  for  burn¬ 
ing  along  a  cylinder).  Therefore,  the  heat  transfer 
rate  near  the  leading  edge  also  increases  as  the  cyl¬ 
inder  radius  decreases.  A  simple  expression  is  de¬ 
veloped  to  evaluate  the  increased  heat  flux,  for  phe¬ 
nomena  2  and  3  above,  by  comparison  to  the  heat 
flux  on  a  flat  plate. 

We  start  with  an  Oseen-type  opposing  flow  having 
a  uniform  velocity  profile  and  constant  properties 
following  development  for  creeping  flame  spread  on 
a  flat  plate  [13,14,17].  The  equation  in  cylindrical 
coordinates  for  a  conserved  scalar,  in  the  gaseous 
flow  is 


df  1  d  / 

dx  ^\dx^  r  dr  \  dr 


(3) 


where  x  and  r  are  the  axial  and  radial  coordinates. 


respectively.  Recall  that  is  a  relative  velocity  into 
the  flame  and  not  a  radial  velocity.  A  methodology 
to  calculate  the  spread  rate  along  the  cylinder  in  a 
manner  similar  to  the  flat  geometry  [18]  is  difficult 
to  implement  because  Bessel  functions  are  needed 
when  taking  the  Fourier  transform,  applying  Wie- 
ner-Hopf  splitting,  and,  finally,  inverting  the  results. 
Instead,  we  evaluate  modification  of  the  heat  trans¬ 
fer  rate  to  the  cylinder  because  of  curvature  by  using 
the  parabolic  form  of  equation  3  away  from  the  lead¬ 
ing  edge  such  that  we  neglect  axial  diffusion,  which 
is  accounted  for  in  equation  1.  Parabolic  equations 
have  been  used  before  in  the  flame  spread  problem 
and  have  been  found  to  give  a  proper  account  of 
energy  redistribution  in  the  gas  and  solid  and  a 
spread  rate  equivalent  to  the  elliptic  problem  [4,17]. 

The  boundary  conditions  are  that  f  has  a  constant 
value  at  the  cylinder  surface,  f  and  is  zero  at 

infinity.  A  closed-form  solution  can  be  found  in  this 
case  [19]  (example  7-18),  which  in  terms  of  the  pres¬ 
ent  formulation  is 


n  h 


]Jcor/R„)Y„{m)  -  ]fw)YJcar/Rj  ^ 

®[/o(®)  +  Y^(®)]  ® 


(4) 


Here,  and  are  Bessel  functions,  with  R,,  the 
cylinder  radius,  and  the  thermal  length,  (5„  is 
Ja^c/V^.  The  derivative  of  equation  4  at  the  cylinder 
surface,  proportional  to  the  heat  flux,  takes  the  form 


d(l> 

dr 


(at  r 


Ro) 


jfi  g-(VHo)W 

(5^  nR„  Jo 


2 

nofffoo)  +  Y^(«)] 

where  standard  properties  of  the  Bessel  functions 
have  been  used. 

An  approximate,  analytical  expression  for  this  in¬ 
tegral  that  is  useful  in  developing  a  closed-form  en¬ 
ergy  balance  for  the  cylinder  similar  to  equation  1 
for  the  flat  surface  can  be  derived  after  finding  its 
asymptotic  expressions  for  small  and  large  3^/R„.  For 
large  radius  (i.e.,  the  flat  plate),  the  integral  equals 
JnR„/2S^  because  (co)  +  Yf  (m)]]  =  1 

for  large  values  of  the  argument,  which  mostly  con¬ 
tribute  to  the  integral.  For  small  radius,  the  integral 
equals 


K 

2 


-  (S^/Rofoj^ 


of  [ln(ta/2)]^ 


which,  following  partial  integration,  is  approximately 
7t/[2  ln(2<5yR„)]. 

An  approximation  of  the  heat  flux  for  any  radius 
can  be  determined  by  combining  the  previous  ap¬ 
proximations  and  checking  with  numerical  integra¬ 
tion  of  equation  5,  which  gives 
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This  expression,  accurate  for  the  Oseen-constant- 
property  flow,  is  useful  in  correlating  experimental 
or  numerical  data.  The  coefficient  C  =  here 
and  the  length  8^  will  depend  on  the  particular 
boundary  layer  flow  along  the  cylinder,  including 
laminar  combustion. 


Heat  Conduction  into  the  Cylinder 

Flame  spread  occurs  along  the  cylinder  because 
the  cylinder  is  preheated  by  the  flame  leading  edge. 
The  sensible  heat  can  be  found  by  solving  the  prob¬ 
lem  for  heat  conduction  into  a  cylinder  for  a  constant 
heat  flux  over  a  specified  preheating  time.  An  ap¬ 
proximation  for  this  preheating  energy  can  be  de¬ 
rived  from  the  exact  solution  in  Ref.  [20]  (p.  194). 
For  the  present  situation,  the  energy  conducted  into 
the  cylinder  for  preheating  takes  the  form 

pMt.  -  r„)VfL,y  (i  - 

Coupling  with  equation  6,  the  flame  spread  equation 
along  a  cylinder  becomes 

pMT.  -  T„)VfL,y  (l  - 


“  — ,  EpL  ErAD  (1) 

^  fi!) 

For  Ljy  greater  than  or  equal  to  R„,  the  cylinder  is 
thermally  thin,  in  which  case  L^y  =  R„. 


Applieations 

For  a  cylinder  at  a  fixed  temperature,  the  numer¬ 
ical  results  for  augmentation  of  the  heat  flux  can  be 
correlated  for  constant  properties  for  several  situa¬ 
tions  [15].  For  a  forced  boundary  layer  flow  along 
the  cylinder  surface,  in  the  parabolic  region  of  the 
flow  removed  from  the  leading  edge  of  the  cylinder 


Nu{R„) 
Mr  (flat) 


3.21 


Rn 


In  1  +  3.21  — 

R. 


with  Nu  the  local  Nusselt  number,  <5^  defined  as  be¬ 
fore,  and  Pr  =  0.7.  For  the  natural  convection  case, 
again  in  the  parabolic  region,  the  constant  3.21  in 


the  expression  for  forced  flow  is  replaced  by  1.80, 
and  the  boundary  layer  thickness  is 

Jg[(TJT„)  -  l]x 

where  T„  is  the  wall  temperature,  and  g  is  the  ac¬ 
celeration  of  gravity.  Thejlifference  in  the  constant 
C  for  the  Oseen  flow  (Jn),  the  forced  flow  (3.21), 
and  the  natural  convection  flow  (1.80)  is  due  to  the 
different  temperature  profiles,  which  also  alter  the 
definition  of  the  thermal  length  From  these  re¬ 
sults,  we  see  that  the  form  of  equation  6  reproduces 
well  the  data  concerning  the  augmentation  of  the 
heat  fluxes  from  the  gas  to  the  solid  because  of  the 
cylindrical  curvature. 


Flame  Spread  on  Cylinders 

The  major  application  of  this  work  is  flame  spread 
along  cylindrical  surfaces.  We  use  numerical  mod¬ 
eling  results  from  both  steady  and  unsteady  com¬ 
putational  codes  [11,12,21].  In  the  numerical  model, 
the  Navier-Stokes  equations  are  solved,  including 
single-step  Arrhenius  gas-phase  combustion  and  sur¬ 
face  pyrolysis.  Surface  re-radiation  and  detailed  gas- 
phase  radiation,  including  feedback  to  the  fuel  sur¬ 
face,  are  accounted  for  by  employing  an  overall 
Planck  mean  absorption  coefficient  in  the  gas  ob¬ 
tained  from  an  overall  energy  balance  consistent 
with  solution  to  the  radiative  transfer  equation  [12]. 
Additionally,  experimental  spread  rate  measure¬ 
ments  from  the  SSCE  for  spread  along  cylindrical 
fuel  samples  in  the  quiescent  microgravity  environ¬ 
ment  of  the  space  shuttle  are  quoted  later  [10,12]. 
The  reader  is  referred  to  Refs.  [11],  [12],  and  [21] 
for  details  concerning  the  numerical  model  and 
properties  employed  and  to  Refs.  [10]  and  [12]  for 
a  detailed  description  of  the  experimental  configu¬ 
ration.  All  results  here  are  for  polymethylmethacry¬ 
late  (PMMA)  as  the  fuel. 

To  apply  equations  1  and  7,  FpL  and  flg  are  needed. 
For  the  flat  surface,  we  notice  that  the  thin  fuel  Epp 
is  independent  of  gas-phase  properties  such  that  it 
can  be  determined  from  solid  properties,  =  1190 
kg/m®,  Cs  =  1465  J/kg  •  K,  =  640  K,  once  Vf  is 
known,  for  which  we  use  results  from  the  numerical 
model.  With  Epp  known,  the  gas-phase  thermal  dif- 
fusivity  that  appears  in  L^y  can  be  determined  from 
the  flat  surface,  thick  fuel  limit,  using,  again,  the 
thick  limit  spread  rate  from  numerical  modeling. 
Obviously,  experimental  values  for  the  Vf  values  may 
be  used  as  well  to  determine  Epp  and  Og. 

For  a  PMMA,  thermally  thin  fuel  thickness  t  of 
0.1  mm,  and  using  equation  1  with  =  T  in  an 
opposing,  forced  flow  of  40  cm/s  at  1  atm,  we  get 
that  Epp  is  127,  305,  404,  and  522  W/m  for  oxygen 
mass  factions  in  nitrogen  of  0.233,  0.533,  0.727,  and 
1.0,  respectively.  The  corresponding  spread  rates  are 
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Fig.  2.  Flame  spread  rate  over  PMMA  flat  surfaces  and 
cylinders  for  =  0.533  in  N2  at  1  atm  and  a  forced, 
opposing  flow  of  40  cm/s  as  a  function  of  thickness  or  ra¬ 
dius.  including  numerical  computations  and  the  analytical 
model  for  cylinders  from  equation  7.  Squares,  numerical 
computations  for  cylinders  from  Ref.  [12];  circles,  numer¬ 
ical  computations  for  flat  surfaces  from  Ref  [12];  solid  line, 
model  equation  for  cylinders. 

2.15,  5.14,  6.81,  and  8.80  mm/s,  respectively  [21]. 
Because  Is  directly  related  to  the  flame  tem¬ 
perature,  the  ratio  of  the  EpL  values  at  different  ox¬ 
ygen  mass  fractions  is  approximately  the  ratio  of  the 
corresponding  oxygen  mass  fractions. 

Corresponding  thick  limit  spread  rates  for  the 
same  oxygen  mass  fractions  are,  in  ascending  mass 
fraction  order,  0.185,  1.18,  1.85,  and  3.42  mm/s. 
These  spread  rates  may  be  used  to  determine  Og 
from  the  thick  formulation  of  equation  1,  and  then 
EpL  and  flg  are  known  for  application  of  equation  7 
at  a  particular  oxygen  concentration,  presuming  C 
can  be  selected. 

The  determination  of  Sg  that  appears  in  L^y  is  com¬ 
plicated  by  the  fact  that  the  correct  velocity  to  use 
forVg  is  an  equivalent  velocity  at  a  distance  Eg  from 
the  surface  that  with  L„  describes  the  velocity  gra¬ 
dient  at  the  surface  rather  than  the  free  stream  ve¬ 
locity  because  of  the  effect  of  the  velocity  profile  on 
distances  over  which  heat  can  be  conducted  [22]. 
The  free  stream  velocity  is  always  higher  than  the 
equivalent  velocity;  cig  then  is  an  effective  property 
such  that  the  ratio  of  flg/Vg  is  identical  to  the  actual 
thermal  diffusivity  divided  by  the  equivalent  velocity. 
Using  the  listed  thick  limit  spread  rates,  and  equa¬ 
tion  1,  we  get  an  of  about  8.80  X  10“^  m^/s  as 
an  average  for  all  the  oxygen  mass  fractions.  The 
equivalent  velocity  for  this  configuration,  from  Ref. 
[22],  is  about  7  cm/s  or  about  6  times  smaller  than 
Vg  so  that  the  actual  thermal  diffusivity  is  around  1.5 
X  10“"*  m^/s.  The  higher  value  is  equivalent  to  a 
thermal  diffusivity  for  1  atm  at  approximately  2400 
K,  while  the  lower  value  is  for  approximately  950  K, 
a  more  realistic  temperature  for  property  evaluation 


[22].  However,  as  long  as  Vg  and  a„  are  used  consis¬ 
tently  in  determining  L^y,  the  result  is  accurate. 

Because  of  the  enhanced  heat  transfer  from  the 
gas  to  the  solid  because  of  the  curvature,  Og  also  ap¬ 
pears  in  the  term  CLg/R,,.  This  a„,  however,  need 
not  be  the  same  as  the  one  used  for  L^y,  which  ap¬ 
pears  in  equation  7  as  a  result  of  in-depth  heating. 
It  is  the  product  Cog  in  the  curvature  effect  that  is 
important  rather  than  Og  alone,  so  that  once  Vg  and 
C  are  selected,  for  use  in  CL^/R„  can  be  deter¬ 
mined  by  requiring  the  thick  limit,  steady  spread 
rate,  if  it  exists,  to  be  recovered. 

Forced  Flow 

In  Fig.  2,  numerical  spread  rate  results  are  pre¬ 
sented  comparing  the  flat  and  cylindrical  geometry 
as  a  function  of  fuel  thickness  and  comparing  the 
numerical  computations  to  the  analytical  model 
from  equation  7,  which  represents  the  behavior  well 
for  the  cylindrical  fuel.  Spread  rates  for  the  cylinder 
are,  as  expected,  higher  than  for  the  flat  surface  as 
the  thickness  decreases  to  the  thin  limit. 

For  Fig.  2,  we  determined  an  Og  of  1.7  X  10“"* 
m^/s  for  use  in  CLg/fi„.  With  the  equivalent  velocity 
being  about  6  times  smaller  than  Vg,  the  actual  Og 
for  use  in  CLg/R,,  would  be  about  30  X  10“®  m^/s, 
which  is  equivalent  to  a  temperature  of  about  350 
K,  which  is  realistic  for  determination  of  heat  con¬ 
duction  into  the  ambient  oxidizer.  Although  the 
value  of  C  used  derives  from  a  boundary  layer  con¬ 
figuration,  the  flame  in  a  forced  flow  is  embedded 
in  the  boundary  layer,  and  the  factor  in  front  of  EpL 
is  to  account  for  the  effect  of  curvature  on  heat  trans¬ 
fer  to  the  surface  rather  than  upstream  heat  transfer, 
which  is  accounted  for  in  the  energy  balance  itself 
through  EpL  and  Eg. 

Figure  3  shows  numerical  computation  and  ana¬ 
lytical  model  equation  results  for  spread  rate  as  a 
function  of  opposing,  forced  flow  in  the  absence  of 
gravity  for  a  fixed  cylinder  diameter  of  3.2  mm,  the 
diameter  used  in  the  SSCE,  using  the  same  model 
equation  deriving  from  equation  7  as  for  Fig.  2.  The 
model  does  not  reproduce  the  computed  spread 
rates  as  well  as  it  does  as  a  function  of  radius  for  a 
fixed  velocity.  The  reason  is  that  once  the  model  pa¬ 
rameters  are  set  for  a  fixed  velocity,  the  behavior  of 
Vf  as  a  function  of  radius  is  accurately  described  by 
the  model  equation,  but  as  the  velocity  changes,  the 
model  becomes  less  accurate  because  the  free 
stream  velocity  should  be  replaced  by  the  equivalent 
velocity  for  accurate  representation  [22].  At  the  low¬ 
est  velocities,  around  5  cm/s,  the  behavior  of  the 
model  and  computed  spread  rates  becomes  dissim¬ 
ilar.  The  reason  for  the  upturn  in  spread  rate  in  the 
model  is  that  C(Lg/R„)/ln(l  +  CL^/R„)  becomes 
large  for  large  LyJR„,  while  E^y]!  —  L,y/2R„)  is 
bounded,  so  a  significant  increase  in  the  curvature 
affecting  the  heat  transfer  from  the  gas  to  the  solid 
is  obtained.  The  computational  results  do  not  show 
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Fig.  3.  Forced  flow,  microgravity  flame  spread  rate  over 
PMMA  cylinders  for  =  0.533  in  Nj  at  1  atm  and  = 
3.2  mm  as  a  function  of  opposing  flow,  including  numerical 
computations  and  the  analytical  model  from  equation  7. 
Squares,  numerical  computations  for  1  atm  from  Ref  [12]; 
solid  line,  model  equation  for  1  atm;  circles,  numerical 
computations  for  2  atm  from  Ref.  [12];  dashed  line,  model 
equation  for  2  atm. 


Fig.  4.  Downward  flame  spread  rates  over  PMMA  cyl¬ 
inders  in  air  at  1  atm  as  a  function  of  radius,  including 
numerical  computations,  the  analytical  model  from  equa¬ 
tion  7,  and  experiment.  Solid  line,  model  equation;  squares, 
numerical  computations  from  Ref.  [12];  circles,  experiment 
[3];  triangles,  experiment  [6];  asterisks,  experiment  [12]. 

the  upturn  because  of  the  effects  of  radiation  de¬ 
pressing  the  heat  transfer  to  the  solid  and  decreasing 
the  gas-phase  thermal  length  scale  [If]. 

Natural  Convection 

Application  of  the  model  equation  in  natural  con¬ 
vection  requires  selection  of  a  velocity  induced  by 
the  spreading  flame,  which  is  [23], 

Using  the  same  properties  as  for  the  forced  flow, 
selecting  the  thermal  conductivity  of  the  gas  fcg  at 
approximately  the  same  temperature  as  the  Og  used 
in  CLg/fi„  in  forced  flow,  and  using  C  =  1.80  from 


natural  convection,  we  get  the  results  shown  in  Fig. 
4  for  spread  rate  as  a  function  of  radius  for  down¬ 
ward  spread  on  cylinders  comparing  numerical 
model  results  with  experiment  and  the  model  equa¬ 
tion.  The  model  equation  describes  the  Vf  versus  R„ 
behavior  correctly,  although  the  Vf  values  tend  to  be 
higher  than  either  numerical  or  experimental  results, 
mostly  due  to  the  approximate  nature  of  the  velocity 
selected. 


Microgravity 

In  quiescent  microgravity,  steady  spread  occurs 
only  for  thin  fuels  [11],  and  so  in  application  of  equa¬ 
tions  1  and  7,  L^y  is  the  fuel  thickness  or  radius.  Ad¬ 
ditionally  in  microgravity,  radiation  from  the  gas  is 
significant  and,  in  fact,  is  the  cause  for  the  occur¬ 
rence  of  only  unsteady  spread  to  extinction  for  thick 
fuels  [11],  and  so  £pL  is  lower  than  the  values  quoted 
above  and  must  be  redetermined.  Fig.  5  shows  nu¬ 
merically  determined  spread  rates  for  flat  surfaces 
at  1  and  2  atm  as  well  as  the  model  equation,  equa¬ 
tion  1,  and  numerically  determined  spread  rates  for 
cylinders  at  2  atm,  along  with  model  results  from 
equation  7  as  well  as  measured  spread  rates  from 
the  SSCE  at  1  atm.  The  computational  data  stop  at 
a  thickness  beyond  which  there  is  no  steady  spread. 
Spread  rates  for  cylinders  are  higher  than  for  flat 
surfaces,  and  the  cylinders  are  able  to  maintain 
steady  spread  for  a  larger  thickness  because  of  the 
enhanced  heat  transfer  that  occurs  as  a  result  of  the 
curvature.  Pressure  effects  for  the  cylinder  appear 
to  be  minimal.  The  behavior  of  the  spread  rate  for 
the  cylinders  is  adequately  described  by  the  model 
equation  with  C  =  3.21  from  the  forced  flow  con¬ 
figuration. 

In  application  of  the  model  equations  to  quiescent 
microgravity,  EpL  was  determined  to  be  164  W/m 
for  the  0.533  mass  fraction  of  oxygen  in  nitrogen  at 
1  atm  (compared  to  305  W/m  in  normal  gravity)  and 
255  W/ m  at  2  atm.  These  values  were  obtained  using 
the  flat  surface  numerically  determined  spread  rates 
for  the  thinnest  fuel  of  Fig.  5,  where  the  spread  rates 
are  highest  and  the  effects  of  surface  radiation  the 
least.  EpL  ™  normal  gravity  is  effectively  indepen¬ 
dent  of  pressure,  while  in  microgravity  it  depends  on 
pressure  because  of  the  decrease  in  total  flame  heat 
loss  by  radiation  with  increasing  pressure.  Although 
absorption  coefficients  increase  with  increasing 
pressure,  the  decrease  in  length  scale  of  the  flame 
results  in  the  reduced  radiation  loss  with  increasing 
pressure  [24].  Epp  also  depends  slightly  on  fuel 
thickness  as  a  result  of  increased  flame  radiation  loss 
with  decreasing  Vf. 

Because  there  is  no  steady  thermally  thick  spread 
in  quiescent  microgravity,  Og  only  appears  in  the 
CLg/fi„  terms.  In  constructing  Fig.  5,  Og  was  set 
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Fig.  5.  Quiescent  microgravity  flame  spread  rates  for  Y„  =  0.533  in  Ng  at  1  atm  as  a  function  of  thickness  or  radius, 
including  numerical  computations,  the  analytical  models  for  flat  surfaces  and  cylinders,  and  experiment.  For  flat  surfaces 
(to  the  left):  closed  circles,  numerical  computations  for  2  atm  from  Ref  [21];  open  circles,  numerical  computations  for 
1  atm  from  Ref  [21];  solid  lines,  model  equation  for  2  atm,  no  radiation;  dashed  line,  model  equation  for  2  atm,  including 
radiation;  solid  line,  model  equation  for  1  atm,  no  radiation;  dashed  line,  model  equation  for  1  atm,  including  radiation. 
For  cylinders  (to  the  right):  solid  line,  model  equation  for  2  atm,  including  radiation;  closed  squares,  numerical  com¬ 
putations  for  2  atm  from  Ref  [12];  dashed  line,  model  equation  for  1  atm,  including  radiation;  X,  numerical  computations 
for  1  atm  from  Ref.  [12];  closed  squares,  experiment  [12]. 


equal  to  a  value  for  ambient  temperature,  which  re¬ 
duces  Lg  compared  to  its  normal  gravity  counterpart 
to  account  for  the  fact  that  radiation  causes  a  de¬ 
pression  in  the  length  scale  such  that  the  mass  and 
thermal  diffusion  scales  are  similar  so  that  steady 
spread  can  be  maintained.  When  the  mass  diffusion 
scale  is  large  compared  to  the  thermal  diffusion 
scale,  steady  spread  is  not  possible  [If], 

An  additional  experiment  executed  in  the  SSCE 
employed  a  hollow  cylinder  of  3.2  mm  radius  and  f 
mm  wall  thickness.  The  model  equation  for  cylin¬ 
drical  spread  gives  a  hollow  cylinder  spread  rate  ap¬ 
proximately  1.54  times  the  solid  cylinder  of  the  same 
radius,  0.71  mm/s  compared  to  0.46  inm/s.  The  ratio 
of  numerically  determined  spread  rates  is  1.36,  0.83 
mm/s  compared  to  0.61  min/s,  and  the  measured 
ratio  from  the  SSCE  is  1.19,  0.93  mm/s  compared 
to  0.78  mm/s.  The  model  and  numerical  results 
oveiyiredict  somewhat  the  curvature  effect  for  the 
hollow  sample.  In  the  model  equation,  gas  curvature 
(the  multiplier  before  £fl)  aiid  radiation  (through 
EpL  and  £rad)  are  accounted  for  separately  when, 
actually,  they  are  not  separate.  Radiation  reduces  Lg 
[24],  which  reduces  the  curvature  effect.  This  cou¬ 
pling  of  gas  curvature  and  radiation  accounted  for 
computationally  but  not  in  the  model  equation  re¬ 
sults  in  overprediction  of  the  curvature  effect  from 
the  model  being  the  largest. 


Conclusions 

Flame  spread  rates  along  cylinders  are  higher  than 
those  over  flat  surfaces  of  the  same  thickness  in  all 
configurations  investigated,  that  is,  forced  and  nat¬ 
ural  convection  flows  and  microgravity.  Analytical 
modeling  shows  that  for  the  cylindrical  geometry, 
there  is  an  enhancement  of  the  heat  transfer  from 
the  flame  to  the  fuel  to  drive  the  spreading  flame 
because  of  the  fuel  curvature,  and  the  rate  of  tem¬ 
perature  rise  of  the  solid  is  faster  because  there  is 
less  material  within  the  cylinder  per  unit  length 
around  the  cylinder  that  needs  to  be  heated  than  for 
the  flat  surface.  The  former  phenomenon  increases 
the  spread  rate  by  a  factor  of 

C(Lg/fiJ/ln(^l  +  C 

above  that  of  the  flat  surface,  and  the  latter  reduces 
the  heated  layer  depth  in  the  cylinder  and  increases 
the  spread  rate  by  a  factor  (1  —  L^J2R„)~^ .  For  the 
thin  cylinder  with  L^y  =  R„,  this  factor  is  two  such 
that  the  cylindrical  spread  rate  is  at  least  twice  that 
for  the  flat  surface.  Erl  is  obtainable  from  numerical 
simulation  or  from  measurement  for  thin  fuels  and 
may  be  thought  of  as  a  material  property,  particularly 
when  flame  radiation  is  not  important. 
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Analytical  expressions  for  cylindrical  spread  de¬ 
rived  from  approximate  solution  to  the  scalar  con¬ 
servation  equation  that  include  the  factors  men¬ 
tioned  above  describe  the  behavior  of  the  spread 
rate  as  a  function  of  cylinder  radius  well  for  forced 
and  natural  flows  and  quiescent  microgravity.  Spread 
rate  variations  with  flow  velocity  are  harder  to  match 
with  the  analytical  expressions  because  of  the  influ¬ 
ence  of  the  characteristic  velocity  gradient  on  the 
appropriate  velocity  to  use  in  the  formulation. 
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